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CO ■ ons in hot hadronic matter using an effective chiral Lagrangian. We find 
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that inelastic processes involving the vector and axial vector mesons reduce 
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, the chemical equilibration time by a factor of ~ 10 compared to the result 



previously calculated in chiral perturbation theory. For a temperature of 
T ~ 150 MeV we obtain a chemical relaxation time of r c /j ~ lOfm/c, which 
is comparable with typical time scales for a hadronic system generated in 



SPS-energy heavy ion collisions. The effect of baryons is also estimated and 
found to be negligible for SPS-energies but important for AGS-energies. We 
predict, that chemical freeze-out should take place at considerable higher tem- 
peratures AT ~ 20 MeV than thermal freeze-out and that the hadronic phase 
would not sustain a pion chemical potential larger than 100 MeV. 
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I. INTRODUCTION 



High energy nucleus-nucleus collisions offer a unique opportunity to explore the large- 
scale properties of quantum chromodynamics, in particular its phase structure at high tem- 
peratures and densities |T| . For sufficiently high collision energy, we expect a phase transition 
into a locally thermalized, deconfined plasma of quarks and gluons, the so called quark-gluon 
plasma (QGP) There may be two different phase transitions, deconfinement and chiral 
phase transition, in QCD at high temperature and/or density. The nature and order of 
the phase transition have been studied as well as possible signatures for the new phase of 
hadronic matter. 

Certainly the system produced in these collisions is not immediately in thermal and 
chemical equilibrium. Secondary interactions among the produced particles are necessary to 
achieve equilibrium. If equilibrium is reached, global observables such as transverse energy 
production can be related to thermodynamic variables, such as energy and entropy density, 
commonly used to characterize these collisions. One can also use the hydrodynamic equa- 
tions to study the evolution of the hadronic system at finite temperature and/or density. 
These conditions considerably facilitate the calculation of certain quark-gluon plasma sig- 
nals, such as the yields of photons, lepton pairs, strangeness and hadrons containing heavy 
quarks. However, even if the heavy ion collision is energetic enough to fully equilibrate the 
system at early time, e.g. in the quark-gluon plasma phase, the system will not necessarily 
remain equilibrated as it expands. Only when the collision rate is greater than the expan- 
sion rate will equilibrium be maintained. It is therefore crucial to understand the following 
questions: How does the system approach equilibrium? What are the time scales for thermal 
and chemical equilibration?, and to which extent is equilibrium maintained in the hadronic 
phase before the system breaks up into the final state hadrons? 

Recently the question of thermodynamic equilibration in the early stages of collisions at 
RHIC and LHC have been studied in the parton model |3|-|5| . In the partonic system created 
early in the collisions, gluons are the dominant components of the system. They reach 
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thermal equilibrium on the time scale 0.3 ~ 0.5 fm/c after collision with an undersaturated 
gluon density. The thermalization of quarks follows that of gluons with the time scale lfm/c 
and reach chemical equilibrium later since the quark pair production by gluon fusion and 
the decay of gluons occur at a somewhat smaller rate. These thermodynamic equilibration 
conditions would be observed in open charm, photon, dilepton, and J /if) production. 

The question of the equilibration in hadronic phase, on the other hand, has been ad- 
dressed by analyzing the observed hadron abundances and spectral distributions f6Hnl . The 
prevailing view requires first a chemical particle freeze-out, from which particle abundances 
are preserved, followed by a thermal freeze-out, after which all interactions between the 
produced particles have ceased. Since the system remains in local thermal equilibrium the 
state can be characterized by a maximum of the entropy consistent with the conservation 
laws for energy, momentum and of the relevant particle number. Neglecting dissipative pro- 
cesses, the corresponding continuity equations then imply that entropy is conserved, i.e. the 
expansion is an adiabatic process. If so, the hadronic system is described by two parame- 
ters, temperature T(x) and a 'chemical potential]]] for each species. Assuming such a 
partial equilibrium, i.e. thermal but not chemical, the observed hadronic abundance can be 
well described. 

As far as strange degrees of freedom are concerned, one expects only a partial saturation 
since the time scale for strangeness saturation is normally larger than the collision time 
scale. The chemical properties of strange particle at freeze-out are regarded as an important 
ingredient to understand the strangeness production in high energy nucleus-nucleus collisions 

E3- 

Naturally, hadronic observables can only tell us something about the conditions at freeze- 
out. In this paper, we are interested in the thermal and chemical equilibrium condition of 



1 A detailed discussion of the meaning of the 'chemical potential' will be given at the beginning of 
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pions in hot hadronic matter at temperatures higher than the freeze-out temperature in order 
to see if the observed freeze-out conditions can be understood within a pure hadronic model. 
We use an effective chiral Lagrangian with vector and axial vector mesons to describe the 
interaction of pions in hadronic matter. In section 2, we first consider the processes that lead 
to thermal equilibrium for pions in hot hadronic matter. Elastic pion collisions, n+n — > n+n, 
turn out to be the principal thermalizing process. Secondly, we clarify the meaning of 
chemical equilibrium and chemical potential of pions in hot matter produced in high energy 
nucleus-nucleus collisions. Thirdly, we estimate the characteristic time scale for chemical 
equilibrium of pions in hadronic matter, first with pions, and then later including resonances. 
Finally, we also explore the effect of baryons on the chemical equilibrium conditions of pions 
in hot hadronic matter. In section 3, the time scales for the chemical equilibrium of pions 
are studied assuming a finite pion chemical potential. 

II. THERMODYNAMIC EQUILIBRATION OF PIONS 

The thermodynamic equilibration in a hadronic system is driven by multiple collisions 
among particles in the system. Generally the equilibration time is directly proportional to 
the collision rate. In an expanding system, therefore, the collision rate should be much larger 
than the expansion rate of the system in order to maintain thermodynamic equilibrium. As 
long as the expansion velocity of the system does not exceed the most probable velocity of 
the thermal particles in the system, this condition is satisfied if the mean free path of a 
particle is shorter than the size of the system. In the following we determine the thermal 
and chemical relaxation time of pions assuming that the system is slightly deviated from 
the equilibrium state, i.e. in the relaxation time approximation. 

A. Thermal equilibrium 

At low temperatures hadronic matter, which mainly consists of pions, can be analyzed 
systematically in the framework of chiral perturbation theory at finite temperature. The 
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low energy theorems, which are obtained from current algebra and partial conservation of 
axial- vector current (PCAC), can be translated into a corresponding set of exact statements 
concerning the coefficients of the low temperature expansion. Since typical pion energies 
are of the order E ~ T at finite T, interactions among pions generate power corrections, 
controlled by the expansion parameter T 2 / f 2 , where f n = 93 MeV is the pion decay constant 
|13| . It is therefore possible to treat the interaction in low temperature pion gas as a 
perturbation. 

The thermal equilibration of pions at low temperatures will be governed by elastic two- 
body collisions, 7r + 7r — ► 7r + n, since inelastic processes contribute in the next leading order 
in the low energy expansion. As long as we consider elastic two-body collisions, the thermal 
relaxation time r t h is given by |14 



J_ _ 1^ _ e -/3£a> 1 f d ^V2 d 3 p 3 d?p A 



T th 4 V J 2E a J (2tt) 3 2£ 2 (2tt) 3 2£ 3 (2vr) 3 2£ 4 

x (27r) 4 5 4 (p a + p 2 - p 3 - p 4 ) 

x £ |A^(7r a 7r 2 ^ 7r 3 vr 4 )| 2 / 2 (l + / 3 )(1 + U), (2.1) 

2,3,4 

where /,'s are Bose-Einstein distribution function. The factor 1/4 in front comes from 
identity of particles at initial and final state and the sum is over the spin and isospin 
degeneracy of particles 2, 3 and 4. The second term (~ e~ f3E ) indicates the contribution from 
inverse reaction. Note that this result is obtained in the relaxation time approximation. In 
the classical limit where Bose-Einstein or Fermi-Dirac quantum statistics is approximated 
by the Boltzmann statistics, this can be written as 

- = ~(1 - e"^) / E <a2 - 34K 2 e"^, (2.2) 
Tth 2 J {Ztt) 6 2 

where a is the cross section corresponding to the relevant reaction and v a 2 is the relativistic 
velocity. 

The mean relaxation time is defined as 



Tth = 77/3 77 \ / d PaTth(Pa)fa(Pa)- (2.3) 

J d Pafa{Pa) J 
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When we use the current algebra result for the ittt scattering cross section we obtain 

12 f 4 

fth « -7^- (2-4) 
With this approximation it has been shown that the mean relaxation time is about 2 fm/c 



at T = 150 MeV [T^J]. Since the mean free path A ~ r in the relativistic limit, the result can 
be compared to the size of the system. Here we assume that the radius of the hot matter 
would be 5 ~ 10 fm. This implies that pions are in thermal equilibrium in hadronic matter, 
mainly due to the two body elastic collisions. 

At temperatures close to the phase transition, other heavy mesons like kaons, vector 
mesons, etc., become increasingly abundant and reduce the thermalization time even further. 
Thus, we conclude that pions can maintain thermal equilibrium in hot hadronic matter 
even at comparatively low temperatures. From formula Eq. ( |2.4| ) we expect a freeze-out 
temperature < 130 MeV for a small system such as S + S or S + Au. This is somewhat 
lower than the values extracted from experiment, which are ~ 150 MeV. The reason for this 
slight discrepancy is most likely due to the flow generated in the reaction, which reduces the 
effective system size. 

We also expect that light vector mesons, p, uj, <p, reach thermal equilibrium in hot 
hadronic matter. For p, uj the dominant reactions will be the collisions with thermal pions 
through heavy resonances such as tt + p — > a\ — > tt + p and tt + uj — > b\ — > tt + uj , see 
Rcf. |15|] . In this reference it has also been shown that the various collisions in hot hadronic 
matter make it possible for phi mesons to be in thermal equilibrium. 

B. Chemical potential of pions in hot hadronic matter 

Before we embark on the discussion of chemical equilibrium we should first clarify what 
we mean when we talk about a 'chemical potential' in the subsequent paragraphs. Strictly 
speaking, a chemical potential is associated with a conserved quantity. In case of pions, the 
conserved quantity would be the charge and, therefore, the only rigorously defined chemical 
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potential would be associated with the charge of the pions (and all other particles involved 
in the system for that matter). This chemical potential, let's call it fi c harge, enters with 
opposite signs in the Bose- Einstein distribution functions for tt + and n~ . In the following 
discussion we will always consider a charge neutral system and, therefore, fj, c harge = 0. 

However, the lifetime of a system created in a relativistic heavy ion collisions is much 
too short in order for the electromagnetic and weak processes to be relevant, so that we 
can safely ignore them from now on. Let us also assume for a moment, that the time scale 
for the strong number changing processes are long compared to the lifetime of the system, 
whereas the elastic processes, which are responsible for the kinetic equilibration, are fast. 
In this case, the number of pions does not change, and, after kinetic equilibrium has been 
established, the pions are distributed according to a Bose-Einstein distribution 

< E ) = e (s-J/r-i ' ( 2 - 5 ) 

where /if is now the chemical potential associated with the conserved pion number. Note, 
that Hf enters with the same sign for all charge states of the pions. In the literature, the 
factor e^/ T is often called a fugacity. Of course in the true thermodynamic limit where the 
system has an infinite amount of time to change the number of particles, /if — 0. Therefore, 
a non vanishing fugacity parameter, \x$ 7^ 0, is often referred to an indication of chemical 
non-equilibrium. This concept can lead to some confusion, as we shall outline now. 

First of all, let us stress the difference between what we shall call true number changing 
processes, such as 

7T7T # tttttttt (2.6) 
and apparent number changing processes such as 

TTTT # p. (2.7) 

As we will demonstrate in the following, the first process drives the fugacity parameter ///to 
zero, leading to an absolute chemical equilibrium. In the second process, on the other hand, 
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the pions are simply 'locked up' in the p-meson, but never really disappear. This process 
ceases to be effective once 2pJ = p^ with pj no necessarily vanishing. Thus, pions and 
p- mesons can be chemically equilibrated with each other at any value of pj. This condition 
we will refer to as relative chemical equilibrium. Consequently, this process and many other 
of that kind, do not drive the system to absolute chemical equilibrium, i.e. pj = 0. 

Whenever in the following discussion we refer to a 'chemical potential' we actually talk 
about the fugacity parameter, p/. Whenever we talk about chemical equilibrium, we mean 
the absolute chemical equilibrium. When particles such as the p-meson and the pions are 
in chemical equilibrium with respect to each other, we will talk about a relative chemical 
equilibrium. After these hopefully clarifying remarks let us now turn to the calculation of 
the chemical relaxation times. 

C. Chemical relaxation time for pions 

The characteristic time associated with chemical equilibrium is determined by a large 
number and variety of inelastic processes. The following processes are subset of possible 
reactions in hadronic matter which change the number of pions; 

(a) 7r + 7T^7r + 7r + 7r + 7r. (2.8) 

(b) IT + 7T # p, 
7T + 7T + 7T # LJ, 

it + p # ai. 

(c) 7T + 7T # p + p, 
7T + 7T # 7T + UJ, 

7r + 7r#7r + ai. (2.10) 

In order to estimate the characteristic time scales for chemical equilibrium we consider 
a system slightly out of equilibrium. For this case the distribution function in phase space 
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(2.9) 



of the various particles takes the form 

where p u is the four-momentum of the particle. X u (x) and /3 u (x) are four vectors constructed 
from the local flow velocity u u (x), the local chemical potential and the inverse local 
temperature f3(x) = 1/T(x): 

\ v (x) = y,{x)u v {x), (2-12) 
(3 u {x) = (3{x)u u {x). (2.13) 

As explained in the beginning, the non-vanishing chemical potential \x is the parameter 
which measures to which extent the system deviates from absolute chemical equilibrium. 

For simplicity we consider hadronic matter in which all thermodynamic variables depend 
only on time and are uniform in space. The change in particle density (n) is then given by 
Boltzmann type rate equations [ |T6|j ; 

^^ = 8n(R mia -R lo8 ^, (2.14) 

where 5n is the number of particles changed in the process and -R ga m, -Rioss are standard 
Boltzmann collision integrals for the production and annihilation rates respectively. 



1. Pion gas 

To be specific, let us consider first hadronic matter at low temperatures which mainly 
consists of pions. Since chiral symmetry suppresses the reaction rates for multi-pion scat- 
tering we consider an inelastic reaction with the minimum number of pions, i.e. the reaction 
7r + 7r^7r + 7r + 7T + 7r. In this case the gain term, i.e. the term that increases the pion 
number is associated with the reaction 7r + 7r^7r + 7r + 7r + 7r and the loss term is related 
with the backwards process. Using the unitarity relation for the matrix element M. 

\M{jXTX — > 7T7T7r7r) | 2 = \M{ixixixix — > 7T7t) | 2 , (2.15) 
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the rate of change for the pion number is given by 

dn n (t) 2 g [ d 3 p 1 f d 3 p 2 f d 3 p 



dt J (2n) 3 2E 1 J (2tt) 3 2E 2 J (2vr) 3 2£ 6 
x (2n) A ^\p 1 +p 2 -p z pe) 

X J2 \M(n 1 H 2 -> 7T 3 - • -7T 6 )| 2 

1,2,. ..,5,6 

x [/i(pi) • • • / 2 (p 2 ) (1 + / 3 (p 3 )) •••(! + /e(P6)) 
- /a(P3) " " " / 6 (p 6 ) (1 + /i(pi))(l + / 2 (p 2 ))] • (2.16) 

Here, the pion density n^if) is the sum over all charge states. The statistical factor S for 
the initial and final state is given by and the sum is over the isospin degeneracy of 
all participating particles. In the Maxwell-Boltzmann limit of the Bose-Einstein quantum 
statistics, we have 

(hu 2(e 2 ^/ T -e 4 ^/ T )/ (T), (2.17) 



dt 



where 



d 3 pi f d 3 p 2 f d 3 p, 



(2tt) 3 2E 1 J (2tt) 3 2E 2 J (2vr) 3 2£ 6 
x (2ir)W 4 \ Pl +p 2 -p 3 pe) 

x J2 |-M(vri7r 2 ^7r3---7r 6 )|V^ + ^/ T . (2.18) 

1,2,. ..,5,6 

Notice, that the number of pions ceases to change if and only if fi n = 0. Thus, the 
process 7r + 7r^7r + 7r + 7r + 7r drives the system to chemical equilibrium and, therefore, 
is a true number changing process according to our previous definition. If we assume that 
the system is only slightly away from chemical equilibrium, i.e. y C 1, we may expand the 
exponentials to first order in Also, in this approximation, for fixed temperature 

n*{t) = nl (l + f^j , (2.19) 

where n° stands for the equilibrium density. The above rate equation Eq. ( ^.17| ) then reduces 
to 
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<^ = -Wo(T). (2.20) 



The solution is given by 



H*(t) = ^(0) exp(-t/r c/l ), (2.21) 

with 

± = (2.22) 

The above approximations are nothing else but the well known relaxation time approxima- 
tion and r c h is called the chemical relaxation time. 

For the process under consideration here the integral Io(T) is given by 

« T > - ™ 6 ), 12 e~, (2.23) 

where o{k\k<i — > 713^4^5^) stands for the isospin averaged cross section and t>i2 for the 
relative velocity. The sum is over isospin of the initial states. In the chiral limit the inelastic 
cross section is |L7 



67 s 3 

J2 ^1*2 -> 7r 3 7r 4 7r 5 7r 6 ) = ^t^^- (2.24) 

With this cross section the chemical relaxation time is proportional to f^/T 9 and estimated 
to be about 26 fm/c at T = 180 MeV (about 40 fm/c with m 7r =138 MeV at the same 
temperature). This is very large compared to the size of the hot matter produced in nucleus- 
nucleus collisions. This implies that the chemical equilibrium in the reaction -kit -Kit-Kit 
would not be reached in the expanding hot hadronic matter consisting of pions only. 

2. Resonance Gas 

In the presence of excited states, i.e. resonances, the number of pions is changed by the 
decay process of resonances (and vice versa), e.g., p # kk, uj # kkk, a% # Kp, and inelastic 
scattering such as pp # kk, kuj # kk, Ka\ # kk etc. 
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(i) Resonance decays 

First we consider resonance decays. In this case the relaxation time will be inversely pro- 
portional to the decay width. Specifically, in case of p tttt we have 

= 2 ( e ^/ T - e 2 ^/ T ) h{p # tttt- T). (2.25) 

at 

Io(p # tttt;T) can be related to the decay width of the p via 



P - (2vr) 3 m p 
« f n p (T), (2.26) 

where the sum is over spin and isospin of vector mesons. f p is the decay width of neutral 
vector meson and n p is the sum over all spin and isospin states. 

Notice, that from Eq. ( |2.25| ), the process p tttt ceases to change the number of pions 
once p p = 2p w . This may very well happen at a finite value p n . Consequently this process 
does not drive the system into absolute chemical equilibrium, as we have defined it. It rather 
ensures a relative chemical equilibrium between p-mesons and pions, which is characterized 
by the condition p p = 2p^. Following our definitions, the process p # tttt is an apparent 
number changing one, since it only shuffles pions back and forth from the p-resonance. The 
fact that the pion chemical potential can take on a finite value means nothing else than the 
conservation of the pion number in this type of process. 

Expanding the distribution functions around the relative chemical equilibrium , one can 
again derive a relaxation time, r^ 1 , which is a measure of how fast the relative chemical 
equilibrium between pions and p-mesons is reached. Following the steps of the previous 
section one finds 

r -=M<TK) =0 ' 4 ~ ' 6fm/c ' (221) 

where n° a indicates the equilibrium density for each species. This is very short, and to a 
good approximation we can assume that pions and p-meson are always in relative chemical 
equilibrium. 
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The same conclusion holds for the decay a± — > irp since r ai ^ p7r = 400 MeV. Thus we 
have p ai = fi n + p p = 3p n where we use fi p = 2p n . However, the decay rate for omega 
mesons into three pions is too small in order to reach relative chemical equilibrium. 

(ii) Inelastic scattering 
Next we consider the inelastic collisions involving resonances such as, tt + tt ^ p + p, 
tt + u) # tt + tt, and 7r + ai # 7T + ii. Again, the rate of change of pion density is given by 
(here in case of tt + tt # p + p) 

dn l = 2 ( e 2 ^ T - e 2 ^ T ) I (n + n^p + p; T). (2.28) 



dt 



As long as the process p ^ tttt is fast compared to the one we are considering here, we can, 
to a good approximation, assume pions and p-mesons are in relative chemical equilibrium, 
i.e. fip = 2p w . With this replacement, the above rate equation Eq. (|2.28|) becomes similar 
to Eq. (|2.17|) and, hence, the chemical relaxation time is given by (see Eq. ( |2.22 )) 

1 4/o(7r + 7r^p + p;T) 



r c h nl 



^2.29) 



The same arguments hold for the process tt + % ^ it + tt, whereas the approximation will 
be not valid in case of7r + u; # tt + tt. For the reactions involving u mesons we assume 
Pu ~ to estimate the corresponding contribution to the chemical relaxation time of pions. 

In order to evaluate the relevant collision integrals, hadronic cross sections have been 
calculated from an effective chiral Lagrangian with vector and axial vector mesons. The 
details are given in appendix A. All processes included in the calculation are summarized 
and results for invariant amplitudes are presented in Appendix B. To take into account the 
off-energy- shell effects we introduce form factors in the calculation of the cross section. The 
standard way to accomplish this is to insert a monopole form factor at each vertex in a 
t— channel diagram 

^ = (2.30) 

where a indicates a species for the exchanged particle. We take a value of A = 1.7 GeV for 
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pjiTi vertices in p exchange reactions |18j and of m p in pion exchange processes. A = m ai 
and 1.7 GeV are taken for aiixp vertices and upn vertices, respectively. 

In Fig. |I| we show the chemical relaxation times for each processes as a function of 
temperature. For comparison the chemical relaxation time corresponding to the reaction 
71717171 7i7i is also shown [17] . The total relaxation time is determined by including all true 
pion number changing processes, with the exception of ti + uj # tt + tt, which is slow because 
of the long lifetime of the u. The result is given by the solid curve. Due to the hadronic 
resonances, the resulting relaxation time is almost an order of magnitude shorter than that 
previously obtained within chiral perturbation theory [[17] . We find that the relaxation time 
is about 1.5 fm/c at T ~ 180 MeV and increases to about 10 fm/c at T = 150 MeV. 

At a temperature of T ~ 150 MeV the total relaxation time is about 10 fm/c, which is 
comparable to typical system sizes created in ultrarelativistic heavy ion collisions involving 
heavy nuclei such as lead. However, the corresponding temperature, for which the thermal 
relaxation time assumes the same value of Tth — 10 fm/c is considerably low, namely about 
110 MeV. From these numbers we expect, that even with the inclusion of the resonance 
chemical freeze-out takes place before the thermal one. Moreover, the above estimate has 
assumed that the resonances are formed instantaneously. Thus, the chemical relaxation 
time, may be slightly larger once the formation of the resonances is taken into account 
properly. 



3. Absorption on baryons 

Next, we take into account the pion number changing processes involving baryons. The 
dominant contributions come from low laying baryons, like nucleon A r (938) and A(1236). 
Since the A resonance decays into Ntt with width T/^ntt = 120 MeV we expect there will be 
a relative chemical equilibrium with respect to the reaction A — > Ntt with //a = A*b + H-n- 
Here \xb indicates the baryon chemical potential, associated with the conservation of the 
baryon number. 



14 



There are also inelastic, i.e true pion number changing, reactions including baryon res- 
onances, such as NN — > iVA. To estimate the characteristic time scale for the baryon 
induced inelastic collisions we use a phenomenological parameterization for the isospin av- 
eraged iVA — > iViV cross section ||19|| ; 

a(NA -> NN) = (pJ/p?)^(VS), (2.31) 

where pf is the momentum in the final NN channel and 

20x( v /i-2.015) 2 
a ^ 3) = 0.015 + 2.015)^ (2 ' 32) 

Here a/s is in GeV, a in mb. Then the time scale T c h is given by 

— = -— e 2 » B/T I B (T), (2.33) 
where we use /ia = Hb + an d 

Here the sum is over the spin and isospin for the initial particles. 

Results are shown in Fig. |[ We use fis = 176, 223, 278, 516MeV for the baryon chemical 
potential which are obtained from the analyses of SPS energy data on strange baryon and 



antibaryon production fqJTT|| and of AGS energy reactions 0. The time scale turns out to 
be about 10 ~ 100 fm/c at T = 180 MeV except for the AGS experiments. This is too large 
to change the number of pions in expanding hot hadronic matter. For AGS experiments 
with fiB= 516 MeV we find that r c h ~ 5 fm/c at T ~ 150 MeV. This implies that the baryon 
effect is important for AGS experiments but negligible in case of SPS experiments. In the 
following we, therefore, ignore the corrections due to the presence of baryons. 

III. CHEMICAL EQUILIBRIUM WITH FINITE PION CHEMICAL POTENTIAL 



So far we have considered the characteristic time scale for chemical equilibration of pions 
when the system is slightly out of chemical equilibrium. This is the case either when the 
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hadronic system evolves from a fully equilibrated quark-gluon plasma phase and hadronizes 
without any changes in chemical properties or when the hadronic system directly is produced 
in thermodynamic equilibrium state. However, if the system is out of chemical equilibrium 
initially, the reaction rate for the pion number changing processes will be affected by the 
finite pion chemical potential. 

We, therefore, extend the previous definition for the chemical relaxation time to the case 
with finite chemical potential and define 

1 ^Sn dn(T, 

T ch ~ [n(T,^)-n°(T)} dt ' l ' j 

which is the ratio of the pion number changing rate over the the excess (or lack) of pions as 
compared to equilibrium. In the limit of vanishing chemical potential, fi n — > 0, the resulting 
chemical relaxation time coincides with that derived in the previous section. 



A. Chemical relaxation time of pions: out-of equilibrium 

In the previous section we found that for a pure pion gas the chemical equilibration rate 
evaluated near chemical equilibrium is too long to to be effective. If the system is far away 
from the chemical equilibrium, then the annihilation (production) reaction rates depend 
on the value of the chemical potential. For the reaction tttttttt — > with finite chemical 
potential of pions, the chemical relaxation time Eq. (|3.1|) is given by 

\ = -^-eWT (1 _ e 2^/T )Io{nm nn . T)/(e M./T _ 1); (3 2) 

Tch n A T ) 

where we take Boltzmann approximation for the distribution of pions. The calculated chem- 
ical relaxation time is shown as dotted line in Fig. |3] for fi w = 100 MeV. When the system is 
out of equilibrium the number changing process becomes considerably faster - about a factor 
of 5 at T = 150 MeV - than in the previous case where the system was near equilibrium 
(note the different scale on the y-axis of Fig. |3|). However, the reaction rate is still large 
compared to the size of the system, r^ h ~ 10 fm/c at T = 180 MeV. Thus, in a hadronic 
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system consisting of only pions, the number of pions ceases to change at the early stage of 
the expansion even if the system is out of equilibrium at the beginning of the evolution. 

Once we include resonances additional processes help to change the pion number. Again 
we assume that the decay processes such as p # tttt and a\ # irp are fast enough to 
maintain the relative chemical equilibrium even if pions have a finite chemical potential. 
From Eq. ( gig) we can see that the pion number is not changed by these reactions. The 



pion number will be changed, however, by the reactions such as pp — > tttt and tt<1\ — > tttt; 
these reactions actually mean that four pions decay into two pions. The chemical relaxation 
time for the corresponding reactions is given by 

-1 = -^r(e WT - e^' T )Upp ^ tttt; T)/(e^ T - 1), (3.3) 

with p p = 2p n . We also get the similar relation for the irai — > tttt; 

1 2 ( e 2 ^ /T -e^ + ^ )/T )I (TTa 1 ^TTTT;T)/(e^ /T -1), (3.4) 



<{T) 
where p ai = 3p n . 

We show the relaxation time with p^ = 100 MeV in Fig. |^. The total chemical relaxation 
time for tttttttt # tttt is given by the sum for all possible channels. With a finite chemical 
potential of pions, the chemical relaxation time of pions now becomes short compared to the 
size of the system. Thus the number of pions would be changed in hot hadronic matter even 
near the thermal freeze-out temperature. The oversaturated pion number will be reduced 
mainly by the reactions pp — > tttt and ttcii — > tttt. Especially axial vector mesons very easily 
interact with pions and annihilate into two pions. However, we should note that the thermal 
freeze-out temperature is also reduced when we include a finite pion chemical potential fli"4 



With non-zero pion chemical potential one also expects an oversaturation of vector and 
axial vector mesons with the relation p p = 2p n and p ai = 3p n , respectively. The above 
result, of course, implies that also the number of vector and axial vector mesons changes as 
quickly as that of the pions. 

In Fig. § we show how the total chemical relaxation time is changed with pion chemical 
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potential at fixed temperature. Even at T = 150 MeV, the pion relaxation time is compara- 
ble to the size of the hot system as long as the pions have a finite chemical potential; r ch as 5 
fm/c at /i n — 50 MeV. If the hadronic system is produced out of equilibrium, for example 
with /i n = 100 MeV, the excess of pions will be reduced by the inelastic reactions involving 
vector and axial vector mesons. These number changing processes will lead to the decrease 
of pion chemical potential and finally cease to be effective as the chemical potential is re- 
duced below a certain value. Therefore, given a freeze-out temperature of about 150 MeV 
we would expect the pion chemical potential not to exceed a value of ~ 50 MeV. 

B. Chemical relaxation time of pions: in equilibrium 

In the previous section we have demonstrated that the chemical relaxation time depends 
on both the temperature as well as the chemical potential. In order to give a realistic 
estimate of the chemical relaxation time, one thus needs to know the chemical potential as a 
function of temperature. The correct treatment, of course, requires the complete solution of 
all kinetic equations including the expansion of the system. This is best done in a transport 
approach and will be addressed in a separate work. In order to give a rough estimate let us 
separate the problem into two pieces. First of all, we determine the chemical potential as 
a function of temperature, assuming the absence of any number changing processes. Given 
that, we can then evaluate the chemical relaxation time using our previous results. 

In order to determine the chemical potential as a function of temperature, we assume 
that the system initially (after hadronization) is in chemical equilibrium. We furthermore 
ignore all true number changing processes. In this case we have only two possible reactions 
in hadronic phase; elastic pion scattering which maintains the thermal equilibrium and the 
decay of the excited states, e.g. p — > im and ai — > irp. The local configuration of pions is 
then described by temperature and a single chemical potential, p. 

Since the system still remains in thermal equilibrium, the state can be characterized by a 
maximum of the entropy consistent with the conservation laws of energy, momentum and of 
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relevant particle numbers. Since we neglect dissipative, i.e. pion number changing, processes 
this implies that entropy as well as the pion number is conserved in the expansion. Thus one 
can assume that the ratio of the effective pion number density n n = n, + 2n p + 3n ai + • ■ ■ to 
entropy density s remains constant during the expansion after hadronization until freeze-out, 

s 

— = constant. (3.5) 
n n 

As a result of keeping s/n n constant, the chemical potential will rise. This is due to the 
overpopulation of pionic states due to the decay of the resonances. It has been estimated 
that pion chemical potential is about 86 MeV at freeze-out temperature [fJ(J]. To be specific, 
we assume that the particles are distributed according to 

fi = T~E> \~~Fn 7' (3-6) 
exp(Ei -m)/T- 1 



where i = ir, p, a x • • • and gi is the spin and isospin degeneracy factor and Ei = yp 2 + mf. 
The entropy density is given by 

_ dP 
S ~df 

[L=const 

where the pressure is 



(37) 



P^/I^M^)}. <-> 

Assuming that the ratio of entropy density of the system to effective pion number density 
is constant throughout expansion, we can obtain the temperature dependence of the pion 



chemical potential as shown in Fig. [5] |20 ] . Here we take two different initial (hadronization) 
temperatures, T h =180 and 200 MeV, where the chemical potential is assumed to be zero, 
fi n = 0. 

With the temperature dependence of the pion chemical potential we can calculate the 
total chemical relaxation time based on the reaction tih-rh — » titt, pp —>■ tttt and a\7i — > tttt. 
We show the result in Fig. |5|. At the initial temperature we, of course, have the same 
value for the chemical relaxation time as that obtained in the limit p n —>■ 0. As the system 
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expands and pions develop the chemical potential the relaxation time becomes shorter than 
that obtained near equilibrium (as shown by dotted line). We find that the r c h is reduced 
by half at T = 150 MeV (r c h ~ 5 fm/c) due to the pion chemical potential. However, even 
with the induced pion chemical potential the chemical relaxation time is considerably larger 
than the thermal relaxation time at the same temperature. 

IV. SUMMARY 

We have studied the thermal and chemical relaxation time scales of pions in hot hadronic 
matter with an effective chiral Lagrangian. From the explicit calculation we show that pions 
in hot hadronic matter are in a phase where elastic collisions rates are very fast compared to 
typical expansion rates of the system. For chemical equilibration the dominant contribution 
comes from the inelastic collision involving ai mesons. Comparing with previous calculations 
fT7|| based on chiral perturbation theory, the inclusion of the resonances has reduced the 
chemical relaxation time by about a factor of 10. When we neglect the formation time of 
these resonances, the resulting chemical relaxation time of pions is 10 fm/c at T = 150 MeV. 
This value is comparable to the size of the hot system produced the collision of large nuclei. 

Given a system size of 5 ~ 10 fm we obtain a thermal freeze-out temperature which is 
small compared to those extracted from experiments |||llj. This might be due to flow effects 
which lead to smaller effective system sizes. If we take the thermal freeze-out temperature 
be about 150 MeV, then the freeze-out size of the system would be 2 ~ 3 fm. On the other 
hand the chemical relaxation time for a system of this size would be T = 170 MeV. This 
implies that chemical freeze-out of pions happens at considerably higher temperatures than 
thermal freeze-out. This result is somewhat at variance with the findings of where 
pion spectra and particle abundances could be reproduced assuming the same freeze-out 
conditions. In order to properly assess the magnitude and importance of this discrepancy, 
a detailed transport calculation including all the number changing processes presented here 
is needed. 
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We also have studied the effect of baryons on the chemical relaxation time of pions. 
Since the effect of baryons is suppressed by their large mass, we consider only low laying 
baryons, iV(938) and A(1236). To estimate the relaxation time we use the phenomenological 
cross section for NN # NA. The effect of baryons is very small and can be neglected in 
SPS experiments. However, it becomes important in AGS experiments where the baryon 
chemical potential is much larger than that in the SPS experiments. 

We have extend the definition for the chemical relaxation time to a system of pions out 
of equilibrium with fi n ^ 0. In this case it is possible that the overpopulated pions will be 
reduced by the inelastic reactions involving vector and axial vector mesons. At T = 150 
with fi n = 100 MeV the relaxation time is about 2 fm/c which is certainly comparable to the 
system size at freeze-out. However, we should note that the thermal freeze-out temperature 
is also reduced when we include a finite pion chemical potential ]14 . 



In order to make contact with reality, we have determined the pion chemical potential 
as a function of temperature assuming isentropic expansion while ignoring number changing 
processes. Based on this relation, we could give a somewhat more realistic estimate of the 
actual chemical relaxation time. We found that in this case the chemical relaxation time 
is about 5 fm/c at T =150 MeV. For systems which are larger than 5fm we, therefore, 
expect that the chemical potential at freeze-out will be considerably smaller than the value 
of 86 MeV, which has been obtained without taking number changing processes into account. 

In conclusion, as long as the effective system size is not considerably smaller than 5 fm, 
a buildup of a pion chemical potential larger than 100 MeV would be very difficult to un- 
derstand. At the same time, we also predict a considerable difference between the chemical 
and thermal freeze-out temperatures. To which extent that is reflected in the data needs to 
be investigated within a transport calculation. 

Future work will concentrate on a transport theoretical calculation of the chemical equi- 
libration. We also plan to extent the present study to include the reactions involving strange 
particles. 



21 



ACKNOWLEDGMENTS 



This work supported by the Director, Office of Energy Research, Office of High Energy 
and Nuclear Physics, Division of Nuclear Physics, Division of Nuclear Sciences, of the U.S. 
Department of Energy under Contract No. DE-AC03-76SF00098. 



22 



V. APPENDIX A: EFFECTIVE CHIRAL LAGRANGIAN 



We consider an effective chiral Lagrangian with vector and axial-vector meson fields 
which are introduced as massive Yang-Mills fields 



C = ^f*Tr[D»UD^} + \f*TrM(U + - 2) 

-^Ti[F^F LtMU + F R F R ^} + m 2 Tr[A%A L >* + A R A R »] 
-iiTi\D^UD u U ] F L ^ + D^D u UF Rtlv } 

+aTr[F^UF R ^U j ], (5.1) 
where U is related to the pseudoscalar fields (f> by 

U = exp 



a/2 ,1 , ^ , r a 



E4>*-7=> ( 5 - 2 ) 

o=l V 2 



U 

and A^(A R ) are left(right)-handed vector fields. The covariant derivative acting on U is 
given by 

D„U = d,U - igAfr - ig\JA% (5.3) 

and Fhj {F R U ) is the field tensor of left(right)-handed vector fields. The A^ and A R can be 
written in terms of vector (V^) and axial-vector fields (A^) as 

A L = l -(V,-A,), A R = ±(V, + A,). (5.4) 



With parameters g = 10.3063, a = 0.3405, £ = 0.4473 |22] we can well describe the properties 
of vector and axial vector mesons. 
From the Lagrangian we have 

Cf vv = |Tr(9 M K - d v VjV»V, 

4 3 ^ = l {y{viTr(d»V„ - d v V^) [A», d v <j>] + ^Tr^A, - d v A,)[d^V v , ^ , (5.5) 

where 5, T]i, and 772 are given by the parameters in the Lagrangian and have the values 
0.347, 0.279, 0.060, respectively. 
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For the four point vertex we have 



4 4 ^ = -^ 2 Tr[y M ,0] 2 



+^Tr[(^K - dyVMPV - d"V^(f> - {dfjYu - d u Vrf<p 2 } 

J TT 

+ ^Tr[(^K - d v V^<f>\y v ,<f>] + [V»AW<P)] 

J TT 

+^Tr\y li ,v v ]& t <f>ar<f>, (5.6) 

J TT 



and also have 



JTT 

T 

da 

J TT 



+ ^TrA^ 



+^Tr[A fi ,d v (f>}[d^,d^ 



+ ^Tr(«Vt - dvA^^U - 4>d^d^), (5.7) 

J TT 

where c's and d's are given by the parameters in the Lagrangian. With parameters given 



in Ref. |22j we have the values 0.00497, —0.00408 for c\ and c 2 , respectively and d\ = 
-2.08, d 2 = -2d!, d 3 = 0.078, andd 4 = -0.0147. 

We also include the gauged Wess-Zumino term in the effective Lagrangian to describe 



the non-Abelian anomaly structure of QCD [23], which leads to an anomalous interaction 



among a pseudoscalar meson and two vector mesons, 

C VVP = -—^—^TT^VudaVpP), (5.8) 

where e^ 01 " is the antisymmetric Levi-Civita tensor with e 0123 = 1. From Eq. ( |5.8| ) we have 

Cunp = -9^ VaP d^ u d aP p ■ 71, (5.9) 

with 



3g 2 

167r2f 



(5.10) 
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VI. APPENDIX B: INELASTIC CROSS SECTION 



In this appendix we summarize the invariant matrix element for each inelastic collision 
considered in section (II. C). The interaction for each process is obtained from the effective 
Lagrangian. 

B.l 7T + TT — > p + p 
B.l.l. 7T+ + 7T- -> p+ + p~; 7T+ + 7T° -> p+ + /9° 



The related diagrams are shown in Fig. 7. The explicit expression for the invariant 
matrix element for each reaction is given by 



M (a) 



A- 



a/2/ (pi+p 2 ) 2 -^ 



1 ^ * 



2m 2 , 



{P1+P2? 



(6.1) 



x 



(Pi - P2) ■ (P3 - P4)fiV + 2 P4 M (P^ - P21/) - 2(pi^ - JtyOPSi/ 



^(psKfa), (6-2) 



[1 + • p 4 - 2^(pi • P4 + P2 ■ P3)]g, 

J TT J TT 



{XV 



-^P^Pzv + 1%{p^P\v + P2»Piv ~ P\»P2v + P2pV\v) U^(p 3 )e u p (p A ), (6.3) 

J TT J TT ) 



where c\ and c 2 are given in previous section. 

The invariant scattering matrix element is then given by the sum as 

M(n + + -R- ^ p + + p~) = M(n + + 7T° -> p+ + p°) 

= .M (cl) + *M (6) + .M (c) . 



(6.4) 



B.l. 2. 7T+ + TT" -> p° + p°; 7T° + 7T° -► p+ + p 



The related diagrams are shown in Figs. 8. The explicit expression for the invariant 
matrix element for each reaction is given by 
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B.3 7T + a\ — > 7T + 7T 



B.3.1 7T+ + a¥ 7T U + 7T + ; 7T+ + 7T+ + 7T~ 



The related diagrams are shown in Fig. 10. The explicit expression for the invariant 
matrix element for each reaction is given by 



M (a) 



JL 



1 



5 



2m 2 p 



(pi - p^y 



X 



Vv^/ A (pi -pa) 2 -m^ V 

(^iP4 • (pi - Pa) - ?72P2 • (Pi - Pa)) (Pi M + Pa M ) 

- ^iP4 • (pi + Pa) - V2P2 ■ (Pi + Pa)) (pm ~ Psm) 



<(P 2 ), 



(6.13) 



.M^ = - 



Mi' 



(c) 



M ic) 



V2J fn (Pi +P2) 2 - mj 



1 - 



2m 2 , 



(Pi+P 2 )' 



x 



V1P1 ■ (Pi + P2) + V2P2 ■ (Pi + P2)) (Pa^ - P4 M ) 

- [vm ■ (Pa - Pa) + i)2ft • (Pa - pS^j (Pa M + P4 M ) 



<(P 2 ), 



1 

7a 

J 7T 



(6.14) 



2d 2 / 2 - 2d 3 p 3 • p A + 2d 4 p 2 • P3^Pi - {d 2 fl + 2rf 3 pi • p 3 + 2d 4 p 2 • P3^P4 
- {^ d ifl ~ d 2fl - 24pi • P4 - 2rf 4 (p 4 - pi) • p 2 )p3 



<(P 2 ), (6.15) 



2 — 



1 

7a 

J 7T 



^2rf 4 / 2 + 2d 3 p 3 ■ p 4 - 2d 4 p 2 ■ p 3 )pi - ^rfx/ 2 - 2d 3 p t • p 3 - 2rf 4 p 2 • p^jp^ 



- (2d 2 fl + Ad 3 pi ■ p 4 + 2d 4 (p4 - pi) ■ P2^P 3 



C(P2), 



(6.16) 



where d^s are given by parameters in the Lagrangian. See Appendix A. 
The invariant matrix element is given by 



7W(tt + + -> 7T+ + 7T") = M (a) + + A< 



(6.17) 
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B.3.2 vr+ + -> 7T° + 7T°; vr+ + af -> tt+ + vr 



r + 



The related diagrams are shown in Fig. 11. The explicit expression for the invariant 
matrix element for each reaction is given by 



M (a) = 



\ 1 



1 



2m 2 p 



(Pi -P3) S 



X 



v 7 ^/ /tt (Pi -P3) 2 -m 2 P \ 
Vm ■ (Pi - Ps) - » • (pi - Pa)) (Pm + Ps M ) 

-(?7iP4 • (Pi +Ps) -» • (Pi +Ps))(Pi M -Ps M ) 



<(P 2 ), (6.18) 



\V2J U{Pl~Pi, 



Pa) 2 — m 2 \ 2m 2 ^ 1 ^ 



x 



1 

J TT 



(^lP3 • (Pi - P4) - W2 • (Pi - P4)) (Pl/x + P4 M ) 

- (»7lP3 • (Pi + P4) - » ' (Pi + P4)) (P1 M - P4 M ) 

2<*i/* - 4/' + 2rf 3 P3 • P4 - 2dt(p 4 + P 3 ) • P^jPt 

- (d 2 fl + 2d 3 p 1 ■ p 4 - 2d 4 pi ■ p 2 )p3 

- (difl + 2rf 3 pi • p 3 - 2rf 4 pi • p 2 )p4 e«! (P2), 

2d 2 / 2 - 4rf 3 p 3 • p 4 + 2rf 4 (P3 + Pa) ■ P^jPi 

- {2difl - 2d sPi ■ P4 + 2d 4 pi • p 2 )p 3 

- (2^ - 2rf 3 Pi • Ps + 2rf 4 Pi • p 2 )p 4 ' < (p 2 ). 



The invariant scattering matrix element is given by 

M(7T + + a' -> 7T° + 7T°) = 7Vf (a) + M (6) + M<i\ 

M(n + + a+ -> 7T+ + 7T+) = 7W (a) + A4 {6) + 7W 2 C) . 



<(P 2 ), (6.19) 



1 

73 

J TT 



(6.20) 



(6.21) 



(6.22) 
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FIGURES 

FIG. 1. Chemical relaxation time of pions in hot hadronic matter. The reactions pp <-> tttt 
(long dashed curve), ttuj <-> tttt (dashed-dot curve), Tra\ <-> 7T7t (dashed curve) and tttttttt <-> 7T7t 
(dotted curve) are considered. The solid curve represents the total relaxation time. These results 
are compared with the thermal relaxation time, r t h, obtained from Eq. (2.4). 

FIG. 2. Contribution of baryon involved reactions to the chemical relaxation time of pions in 
hot hadronic matter. Results are obtained with ps = 176 MeV (dotted line : SPS), 223 MeV 
(dashed line : SPS), 278 MeV (long-dashed line : SPS) and p B =516 MeV (solid line : AGS). 

FIG. 3. Chemical relaxation time of pions in hot hadronic matter with finite pion chemical 
potential. The contributions from the reaction pp <-> tttt (long dashed curve), ttcl\ <-> tttt (dashed 
curve) and tttttttt <-> tttt (dotted curve) are shown. The solid curve is for the total relaxation time. 
The results are obtained with p w = 100 MeV. 

FIG. 4. Chemical relaxation time of pions in hot hadronic matter as function of pion chemical 
potential. We consider four different temperatures; from top T = 140, 150, 160, 170 MeV. 

FIG. 5. Temperature dependence of the pion chemical potential: We assume two different 
initial (hadronization) temperatures, = 180 (lower curve) and 200 MeV (upper curve). 

FIG. 6. Chemical relaxation time of pions with pion chemical potential shown in Fig. 3. Here 
we assume that 1)^=180 MeV. 

FIG. 7. Diagrams for the reaction tt + + tt~ — > p + + p~ and tt + + tt° — > p + + p°. Dashed line 
and solid line indicate the pion and p meson, respectively. 

FIG. 8. Same as Fig. 7 for the reaction tt + + tt~ — > p° + p° and tt° + tt° — > p + + p~ . 

FIG. 9. Diagrams for the reaction ttuo — > tttt. Dashed line, solid line and double solid line 
indicate the pion, p meson and to meson, respectively. 
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FIG. 10. Diagrams for the reaction tt + + a® — > ir° + 7r + and 7r + + a 1 — ► 7r + + 7r Dashed line, 
solid line and double solid line indicate the pion, p meson and a± meson, respectively. 

FIG. 11. Same as Fig. 10 for the reaction ir + + — ► 7r° + 7r° and 7r + + — ► 7r + + ir + . 
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